Evidence of a subenergy gap in the overdoped regime of Yi_ x CaL x Ba2Cu 3 7 s thin 

films from THz Spectroscopy 



N. Bachar, 1 * 2 ^ E. Farber, 1 E. Zhukova, 3 B. Gorshunov, 3 and M. Roth 4 

1 Laboratory for Superconductivity and Optical Spectroscopy, 
Ariel University Center of Samaria, Ariel, 40700, Israel 
2 The Raymond and Beverly Sackler School of Physics & Astronomy, Tel Aviv University, Tel Aviv, 69978, Israel 
3 Submillimeter Spectroscopy Department, Prokhorov General Physics Institute, Russian Academy of Sciences, Moscow, Russia 

4 Department of Applied Physics, Hebrew University, Jerusalem, Israel 
(Dated: March 4, 2013) 

We measured the terahertz (THz) complex conductivity ol Ca doped YBa2Cu307_<5 thin films 
in the frequency range of 0.1 to 3 THz (3 to 100 cm" 1 ) and at a temperature range of 20 to 
300 K. The films were measured using both time domain and frequency domain THz methods. We 
showed evidence for the existence of a sub-gap in overdoped Yi_ a; Ca a ;Ba2Cu307_<s samples doped 
with 5% and 10% Ca. Evidence for the opening of this sub-gap appears as a sharp decrease in the 
spectrum of the real part of conductivity at frequencies equivalent to a gap energy of 1 meV and 
is more prominent with increased doping. This decrease in conductivity can be explained by using 
d x 2_ y 2-wicve pairing symmetry with an imaginary part of is or id xy which suggests node removal. 



Symmetry of the order parameter in high-T c cuprates 
superconductors and the existence of nodal lines along 
the diagonal of the Brillouin zone, the TX direction, has 
been and is still one of the theoretical and experimental 
questions to be addressed 1 -. There is a degree of con- 
sensus that the pairing symmetry consists of a dominant 
d x 2_ y 2-wave where there is an extra imaginary compo- 
nent which is doping and magnetic dependent. This ex- 
tra component is still puzzling and leads us to important 
questions that remain to be resolved. One of these ques- 
tions is whether the additional imaginary component of 
the order parameter (OP), associated with surface cur- 
rents^, occurs only at the surface or whether it is in fact 
a bulk property 3 -. While a surface-associated imaginary 
component would still be compatible with a pure d x 2_ y 2- 
wave bulk OP, a bulk- associated imaginary component 
would not be. This distinction has important implica- 
tions for our understanding of the high T c superconduc- 
tivity mechanism. Some of the proposed models require 
pure bulk d-wave symmetry 4 ^—, while others do noi£~— . 

Unlike conventional superconductors, in high- cuprates 
there is much research which supports anisotropic sym- 
metry exhibiting nodes in the order parameter, i.e. 
d a; 2_j / 2-wave symmetry, see for instanc e 11 ' 12 . In this case, 
the d a; 2_j / 2-wave superconducting gap is estimated to be 



A fa 20 meV for YBa 2 Cu 3 07_ 



Nevertheless, possi- 



ble mixed states of complex order parameter d x 2_ y 2 + is 
or d x 2_ y 2 + id xy are also allowed within the symmetry 
group of YBa 2 Cu307_5 — . Tunneling^ - — and penetra- 
tion depth-?- measurements show an additional imaginary 
component estimated to be of the order of 8 ~ 2 meV. 
This component was induced by magnetic field in opti- 
mally doped YBa2Cu307_«5 films, and observed in over- 
doped Yi_ x Ca x Ba2Cu307_,5 films without applying an 
external magnetic field. 

Infrared spectroscopy, which was used extensively in 
the field of high-T c cuprates, holds the same energy scale 
for the superconducting energy gap and the excitations 



energies which are assumed to be relevant for the forma- 
tion of superconductivity^. However, by using Far In- 
frared techniques such as THz spectroscopy, one is able 
to obtain the quasiparticles's dynamics far from these ex- 
citations energies spectrum. The vast use of THz spec- 
troscopy in the past years has given us the advantage 
of probing bulk material properties while achieving a 
wide frequency range. By these means, we can focus on 
probing additional gap contributions in the THz range 
where the energy scale is of the order of 1 meV. Our 
work is focused on probing the complex conductivity of 
Yi_ a; Ca a ;Ba2Cu307_5 in the THz frequency range where 
the energy scale of the radiation is equivalent to the en- 
ergy scale of the sub-gap, 8 ps 2 meV. As will be further 
elaborated, we tried to estimate a unique frequency and 
temperature dependence of the complex conductivity in- 
dicating the order parameter symmetry in our samples. 



I. EXPERIMENTAL 

In the current research two batches of 
Yi_ a; Ca a ;Ba2Cu307_«5 thin films were deposited. 
The films were grown to a thickness of 500 - 600 A on 
LaA103 substrates using 10% and 5% Ca. The films 
were deposited by a DC sputtering method at a growth 
temperature of 820 oC, using a pressure template proce- 
dure^ including additional water vapor of 2 - 3 mTorr. 
XRD measurements show a well oriented c-axis films. 
The samples' transition temperature was characterized 
using resistance and induction measurements. The 10% 
films showed typical T c of 78 K ± 1 K while the 5% Ca 
doped films showed typical T c of 84 K ± 0.7 K. From the 
curvature of the resistivity as a function of temperature, 
we have deduced that the samples are in the overdoped 
regime. This is done by taking the second derivative of 
the resistivity curve as a function of temperature. 

The first batch of films was measured using a Time 



2 



Domain THz Spectroscopy method applying a TeraView 
TPS Spectra 1000™system. Time Domain THz spec- 
troscopy method has been reviewed in several paper o 20 ' 21 . 
The complex transmission spectra of the films, substrate 
and reference (clear aperture) were recorded from room 
temperature down to 20 K with steps of about 1 K and 
in the frequency range of 0.1 - 3 THz (3 - 100 cm -1 ). 

The second batch of films was measured by using a 
THz Quasi-Optical method applying an Epsilon coherent 
source spectrometer. This is a well-known quasi optical 
method, which has been reviewed and used by several 
groups 2 ^. The absolute transmission function and the 
transmission phase shift of the films, substrate and ref- 
erence were recorded from room temperature down to 
20 K in steps of 50 K above 100 K and steps of 10 K 
below 100 K. The data was recorded for four Backward 
Wave Oscillators THz sources in the frequency ranges of 
2 - 4 cm" 1 (QS3-120 source), 8 - 23 cm" 1 (QS1-370 and 
QS1-700 sources), and 28 - 41 cm" 1 (QS1-1100 source). 

Analysis of the complex transmission in both meth- 
ods was first performed for the bare substrate. The sub- 
strate exhibited unfamiliar sharp absorption at frequency 
of about 65 cm -1 due to interaction of polarization with 
the a-b twins which were not oriented with the sample 
edges. Therefore, all samples, i.e. films and substrate in 
both methods, were rotated by approximately 45 degrees 
to the polarized beam, in order to reduce this effect. In 
this case, we were able to measure our samples with no 
absorption in the entire frequency range. It is important 
to mention that the absorption did not affect the trans- 
mission at lower frequency range below 40 cm -1 . The 
data for the LaA103 substrate was fitted using Lorentz 
oscillator for the phonon absorption with the resonance 
frequency of 190 cm -1 — and resulted in an almost con- 
stant dielectric function of E\ « 24 in the frequency range 
1 to 60 cm -1 at temperatures 20 to 300 K. 

For data analysis of our superconducting film, we used 
the two-fluid model. The scattering rate and plasma fre- 
quency of the quasi-particles and plasma frequency of 
the superconducting carriers were determined by a least- 
square fitting of the transmission and phase spectrum 
from which the complex conductivity was deduced 2 ^. In 
order to achieve the best fit, we have minimized the good- 
ness of fit parameter, % 2 , by additionally tuning our fit 
parameters as a function of frequency and applying the 
Variable Dielectric Function 2 ^ for the THz-TDS data and 
Fresnel equations^ for the Quasi Optical system data. 



II. RESULTS AND DISCUSSION 

Two batches of 5% and 10% 
Yi_ 2; Ca 2 ;Ba2Cu307_5 thin films were measured utilizing 
both of the spectroscopy methods. The transmission of 
the films decreases with temperature as shown in Fig- 
ure [1] and in Figure 2(a) for time domain and frequency 
domain measurements, respectively. Below T c 
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FIG. 1. Typical THz pulse measured in transmission through 
Yi_ a; Ca I Ba2Cu307_5 (x=10%) thin film at temperatures of 
300 K and 25 K. Pulse exhibit reduction in amplitude and 
shift in position, i.e. time. Below T c « 78 K The com- 
plex transmission function is obtained by the ratio of pulse 
transmitted through sample to a pulse transmitted through 
reference. 



is also a change in transmission phase shift which is 
shown in Figure [Tj and in Figure 2(b) In the case of the 



frequency domain measurements, these drastic phase 
shift changes are observed mainly in the low frequency 
regime. 

We have used the transmission function of four inter- 
faces, i.e. sample with two layers. For the data analysis, 
the two-fluid model was used in order to describe the 
complex conductivity of the thin film 



<7i = <7 n (oj,T) + <t s (w,T) 



(1) 



Where a n is the complex conductivity of the normal car- 
riers and a s is the complex conductivity of the supercon- 
ducting carriers which are described as 



<7„(w,T) = oj Pn (T)[ 



T(T) 



r 2 (T) + w 2 r 2 (r) + , 



& s (uj,T) = w Ps (T)[6(uj = 0)+i- 



(2) 



(3) 



along 

with a sharp decrease of transmission amplitude, there 



where T(T) is the normal carriers' (or quasiparticles' at 
temperatures below T c ) temperature dependent scatter- 
ing rate, 0J Pn and uj Ps are the plasma frequency of normal 
and superconducting carriers, respectively. 

The data analysis was confined to a change in the scat- 
tering rate of the normal carriers or quasiparticles below 
T c and change in the plasma frequencies for both types 
of carriers. Our analysis of the complex transmission 
data results in the fitting parameters' temperature de- 
pendence as shown in Figure [3] In the following section 
temperature and frequency dependence of the complex 
conductivity will be analyzed applying the frequency in- 
dependent parameters of the two-fluid model followed by 
a frequency dependent fit. Finally, the anomaly in the 
complex conductivity seen at low frequencies and well 
below T r will be reviewed. 
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FIG. 2. The Transmission (a) and Transmission Phase Shift 
(b) of a Yi_ 2! Ca a; Ba2Cu307_5 (x=5%) thin film measured 
with the Quasi Optic system for selected temperatures and 
frequency ranges. Figure (a) shows decrease of transmission 
as cooling down and the decrease of transmission particularly 
at low frequencies below the critical temperature. Figure (b) 
shows change in phase shift amplitude and increased oscilla- 
tions at low frequencies and below the critical temperature. 
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(b)Plasma frequency of the superconducting 
carriers 

FIG. 3. Frequency independent superconducting carriers 
plasma frequency and quasi-particles scattering rate obtained 
by fitting to the two fluid model using a delta peak for the Su- 
perconducting (SC) carriers and Drude contribution for the 
quasi-particles. The figure shows a drastic decrease in the 
scattering rate at T c while saturating at low temperatures. 
At T c there is also an increase in oj Ps due to the formation of 
pairs. 



Above T c the temperature dependence of complex 
transmission is denned mainly by a decrease in its ab- 
solute value. This behavior can be understood by a de- 
crease in the normal carriers' scattering rate in the two- 
fluid model as shown in Figure 3(a)| As a result, the 
real part of conductivity, ai(u, T) , slightly increases as 
temperature decreases. This increase of a\(uj,T) above 
T c is not followed by any change in the imaginary part 
of conductivity, 02 (w, T) . 

Below T c there is a sharp decrease of transmission and 
in addition, a drastic change in the phase shift. The sharp 
decrease in transmission and the change in phase shift are 
observed mainly at low frequencies. This behavior is well 
explained by the formation of superconducting carriers 
at the transition. Therefore it is consistent with the in- 
crease in the superconducting plasma frequency as shown 
in Figure [3(b) | In addition, the scattering rate decreases 



sharply as shown in Figure 3(a) This change is also 



related to condensation of normal carriers into the su- 
perconducting state which results in a weakening of scat- 
tering mechanisms in the film. As a result of the above 
changes, our analysis shows a sharp increase in both the 
real and imaginary parts of the conductivity as shown in 



Figure [4] and Figure [5] for the frequency domain and time 
domain methods, respectively. It is important to note 
that the superconducting contribution is related mainly 
to the imaginary part of the conductivity (Eq. [3]) in our 
frequency range, therefore explaining the drastic change 
in phase shift at low frequencies. According to Eq. [3l the 
imaginary part of the conductivity, 02(0;, T) , in the two- 
fluid model is diver ging a t low freque ncies; this tendency 
is shown in Figure |4(b)| and Figure 5(b) In order to 
understand this behavior, the change of the conductivity 
spectral weight below T c has to be considered. In the nor- 
mal state and due to a high scattering rate, the real part 
of conductivity shows an almost frequency independent 
behavior in our frequency range. However, below T c , the 
scattering rate T decreases drastically therefore quench- 
ing the spectral weight down to lower frequencies. The 
change in the spectral weight at the transition is shown 
in our data below T c where a\{oj,T) exhibit a gradual 
change as a function of frequency. 

At lower temperatures below T c , our analysis shows 
that the scattering rate is maintaining somewhat fixed 
value, consistent with other reports^. In this tempera- 
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(b)Imaginary part of conductivity - a2(u>,T) 

FIG. 4. The real part of conductivity as a function of 
frequency (b) The imaginary part of the complex conduc- 
tivity as a function of frequency. Both figures represent 
Yi_ a; Ca a: Ba2Cu307_i (x — 10%) thin film measured by the 
Quasi Optic method. Figure (a) shows selected temperatures 
above and below T c . Below T c the real part of conductivity 
shows a non monotonic behavior with a decrease at low fre- 
quencies. Figure (b) shows the increase of the imaginary part 
as 1/w according to Eq. [3] 



ture range and due to additional spectral weight removal 
from the Drude part of the quasi-particles to the delta 
peak of the superconducting carriers at zero frequency, 
the plasma frequency of the quasi-particles also starts 
to decrease. This behavior is the main cause for the 
"coherence-peak" -like behavior observed in the real part 
of the conductivity for optimally doped YBCO films 28 
and in our current research for overdoped YBCO films 29 . 

After applying the frequency independent parameters 
of the two-fluid model as a first approach, we can deploy 
other techniques in order to tune our results for achiev- 
ing the best fit. In the frequency domain we have just 
used the Fresnel' equations 2 ^ for complex transmission of 
a two layer sample and in the time domain we have used 
the Variable Dielectric Function 2 ^ for the film's data. In 
the frequency domain, the complex conductivity was cal- 
culated by applying the Fresnel equations on each data 
point without using Kramers Kroning relations. In the 
time domain, the VDF function is modeled as a sum of 
Lorentz-Drude oscillators fixed at each data point which 
were tuned in order to achieve the best fit. These ap- 




(a)Real part of conductivity - u\ (cj, T) 




(b)Imaginary part of conductivity - <J2(u>, T) 



FIG. 5. The real and imaginary part of the complex conduc- 
tivity for Yi_ a; Ca a ;Ba2Cu307_j (x — 10%) thin film mea- 
sured by Time Domain Spectroscopy system. Figure (a) 
shows onset frequency in which the real part of the conduc- 
tivity increases below T c , and followed by a decrease at about 
10 cm -1 and at temperatures below 40 K. Figure (b) shows 
the increase of the imaginary part having 1/uj dependency 
according to Eq. [3] 



proaches, in both techniques, enable us to obtain bet- 
ter agreement between our data points and a customized 
fit model. It is important to note that the parameters 
r,u> Pn ,ujp n which were used in the first fit, and shown in 
Figure |3l were defined as frequency independent param- 
eters. As a result, for frequencies much smaller than T, 
the real part of conductivity should be almost frequency 
independent. However, after applying the VDF or solv- 
ing the Fresnel equations, one obtains a complex con- 
ductivity which deviates from the frequency-independent 
model. 

Following the above-mentioned procedure, further 
analysis of the complex conductivity data below T c shows 
some known features for high-T c thin films. While 
02(0;, T) gradually increases due to the increase in the 
superconducting carriers' concentration as temperature 
drops, <j 1 (u>,T) , at a fixed frequency, reaches a maximum 
at about 40 K. This is a well-known signature for high 
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FIG. 6. Imaginary part of conductivity multiplied by the 
frequency for 10% Ca doped sample. The imaginary part 
shows a deep at frequencies equivalent to the onset in the real 
part 



FIG. 7. Simulation of the DOS as function of energy for pure 
d x 2_ y 2-wave and complex energy gap symmetries. The figure 
shows a cutoff energy scale, S, underneath no quasiparticle 
states are available 



quality thin films. The maximum conductivity as a func- 
tion of temperature can also be related to the competition 
between changes in the scattering rate and the plasma 
frequency of the quasi-particles as temperature reduces. 
As a function of frequency, the real part of conductivity 
should increase or at least retain a maximum value at 
low frequencies. However, our data show some anoma- 
lies in the complex conductivity at low frequencies which 
was not reported for optimally doped YBa 2 Cu307_,j . 
Following this temperature and frequency dependence of 
ci (u, T) , there is a sharp decrease of a\ (w, T) at low fre- 
quencies, which becomes stronger as temperature drops 
down to 20 K. This decrease is shown in Figure 4(a) and 
Figure 5(a) at a frequency of about 10 cm" 1 . The on- 



set drop in a\(uj,T) is observed as a dip in wo^w, T)for 
these overdoped samples at about the same frequency as 
shown in Figure [6] As shown in Eq. [3j the dominating 
part in the imaginary uht2(u), T)is proportional to the su- 
perconducting carriers plasma frequency. 

The uncommon feature observed in the real part of 
conductivity is inconsistent with a pure d x i _ y i -wave su- 
perconductor analyzed by the two fluid model. Instead 
of showing a constant conductivity value at low frequen- 
cies, as expected from the Drude function, the real part 
exhibits a drop below 10 cm" 1 , as shown for the 10% 
Ca overdoped samples. This feature appears at temper- 
atures below 40 K and seems to be enhanced as temper- 
ature decreases. 

Comparing our data to previously reported data in 
the same frequency range using optimally doped sam- 
ple d 30 ! 31 one can see that this non-monotonic behavior 
of cti(w,T) at low frequencies and below T c is reported 
here for the first time. This behavior is observed for 
the first time in our current research on Ca overdoped 
YBa2Cu307_<5 thin films. 

The current paragraph will elaborate on some possible 
explanations of our data. We can regard this non mono- 
tonic behavior as result of a frequency dependent scat- 
tering rate or plasma frequency which deviates, mainly 
at low frequencies, from the values shown in Figure [3] 
A decrease in the normal part of conductivity can be 



explained as an increase in the scattering rate or as a de- 
crease in the quasiparticles plasma frequency at the fre- 
quency range of 10 cm -1 . The increase in the scattering 
can be compatible with Ca impurities causing localiza- 
tion at low frequencies. However, it is not at all clear 
why it should appear only at temperatures quite below 
T c . Instead, we suggest that this decrease in conduc- 
tivity should be considered as a decrease in the density 
of states (DOS). The decrease in the plasma frequency 
at low frequencies is just a direct observation of this be- 
havior. In the case of YBa2Cu307_,5 we have compared 
the pure d x 2_ y 2-w&ve DOS to pairing symmetries such as 
d x 2_ y 2 + is and d x 2_ y 2 + id xy , i.e. mixed d x 2_ y 2-w&ve . 
Figure [7] shows that for the pure d x 2_ y 2-w&ve there are 
states down to zero energy. However, for the complex 
order parameter there are no states available below a 
sub-gap. We suggest that this sub-gap is a result of node 
removal in the pure d x 2_ y 2-w&ve by an additional imag- 
inary d xy or s component. The sharp decrease in the 
density of states should be observed as a decrease in the 
quasi-particles concentration or, experimentally, by a de- 
crease in their plasma frequency As the real part of con- 
ductivity is a direct probe of the quasi-particles' plasma 
frequency, a decrease in its value is observed as this sub- 
gap in energy. 

The idea of an additional imaginary part in a super- 
conducting order parameter was used for calculating the 
complex conductivity by Schiirrer et al? 2 . In their work, 
s + id pairing symmetry was assumed using the Eliash- 
berg equations. In their report, several quasi-particle 
density of states (QDOS) spectra were calculated for dif- 
ferent imaginary and real part ratios. In the case of 
small s compared to the d part and by assuming that 
the QDOS yields the absolute value of the order param- 
eter, the results of Schiirrer et al.&re consistent with our 
data. The complex conductivity which was calculated 
by Schiirrer et aLusing a s + id pairing symmetry shows 
several striking features. The real part of conductivity, 
ui(uj,T) , shows an onset followed by a decrease at low 
frequencies which are equivalent to the sub-gap energy 
scale observed in our current work. In addition, the imag- 
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inary part product, (T2 (w,T) , shows a dip developing at 
the same energy scale as the onset frequency in a\(u), T) . 
The complex conductivity data in our films, which shows 
deviations from the pure d x 2_ y 2-wave and the frequency 
independent parameters of the two fluid model analysis, 
exhibits the same features as described above. We believe 
that the only possible interpretation of these deviations 
can be achieved by an additional imaginary part of the 
order parameter, as suggested by Schiirrer et al. 

Finally, we would like to clarify at this point that the 
reduction of cti(cj,T) at about 10 cm -1 , which is as- 
sumed to be the sub-gap frequency, is at a much lower 
energy than the imaginary sub-gap measured directly by 
tunneling or obtained indirectly by microwave London 
penetration depth to have a value of about S » 2 meV(~ 
16 cm _1 )ii. We cannot rule out a possible explanation 
for this discrepancy where this extra imaginary gap com- 
ponent is a surface effect and reduces once bulk is exam- 
ined. 



III. CONCLUSIONS 

The complex conductivity of Ca doped 
YBa2Cu307_5 thin films was measured by both 
frequency and time domain spectroscopy techniques 
in the THz frequency range. We have shown distinct 
evidence for a sub energy gap in the THz frequency 
range which develops below T c and exhibits itself as a 
significant decrease in the real part of conductivity with 
a typical onset frequency of 10 cm -1 . The sub-gap 



in the real part of the conductivity shows up as a 
dip in aj<7 2 (w,T)at the same frequency, representing 
a decrease in the superconducting plasma frequency. 
These features in the measured complex conductivity 
spectrum match the calculated spectrum by Schiirrer et 
abusing a complex order parameter with imaginary part 
of 25 ss 1.2 meV . We suggest that this gap is related to 
node removal in the overdoped regime of YBa2Cu307_a . 
This idea has previously been reported utilizing different 
surface sensitive tunneling experiments and microwave 
indirect measurements. However, this gap appears at 
lower energies than in previous research. At the current 
stage we cannot rule out the effect of Ca atoms, which 
might act as impurity centers causing localization and 
this could be done by a future research in which the 
Ca doped samples are annealed into the optimal doped 
regime. In order to avoid the influence of a-b twins on 
frequency response, caused by the substrate, which is 
presumably temperature independent, overdoped thin 
films, should be grown on un-twinned substrates. 
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